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Abstract 

We  revisit  the  problem  of  evaluating  the  performances  of  communication  systems  with 
optical  amplifiers  and  a  wideband  optical  filter.  We  compute  the  exact  probability  of  error  and 
the  optimal  threshold  and  compare  them  with  those  predicted  by  Gaussian  approximations  for 
ASK,  FSK  or  DPSK  modrdations,  both  for  ideal  photodetectors  and  for  the  case  where  shot 
noise  is  significant. 

1  Introduction 


We  consider  the  problem  of  receiving  modulated  signals  that  have  passed  through  optical  ampli¬ 
fiers.  An  optical  amplifier  with  power  gain  G  can  be  modeled  as  a  linear  optical  field  amplifier 
together  with  a  source  of  white  Gaussian  noise  over  the  bandwidth  of  interest.  The  noise  has  a 
two  sided  spectral  density  Ao/2,  with  Nq  =  Ngphi/{G  —  1)  [21].  There  h  is  Planck’s  constant,  i/ 
is  the  frequency  of  interest,  and  N,p  is  a  factor  that  takes  various  imperfections  into  account, 
its  value  is  ideally  1. 

In  the  case  of  amplitude  shift  keying  ( ASK )  the  signal  is  either  0  or  an  optical  pulse  of  dura¬ 
tion  T.  A  possible  receiver  structure  consists  of  a  jmlarizer  passing  only  the  signal  component, 

an  optical  bandpass  filter  of  bandwidth  B,  a  photodetector  and  an  integrator  over  time  T. 

*'riie  authors  gratefully  aekii  owl  edge  suiiport  for  this  researeh  from  NSF  (grant  SS(I2991-N<  T ).  and  DARPA 
(grant  i’l<.)(52S-90-(:'-0(,)02) 
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In  the  case  of  frequency  shift  keying  (FSK)  there  are  two  orthogonal  signals,  two  filters  whose 
frequency  response  (each  of  bandwidth  B)  do  not  overlap,  two  photodetectors  and  a  circuit 
integrating  the  difference  of  the  photodetector  currents  over  time  T.  Very  similar  structures 
can  be  used  for  other  binary  orthogonal  signals,  such  as  in  polarization  modulation  and  in 
binary  pulse  position  modulation. 

In  the  case  of  differential  phase  shift  keying  (DPSK)  the  signal  during  a  bit  time  is  a  replica 
of  the  signal  sent  during  the  previous  bit,  with  phase  possibly  shifted  by  t.  The  receiver  consists 
of  an  optical  filter  of  bandwidth  i7,  a  Mach-Zehnder  filter  with  differential  delay  T  (adjusted 
so  that  vT  is  an  integer),  two  photodetectors  (one  for  each  output  branch  of  the  Mach-Zehnder 
filter)  and  the  same  photodetectors  and  integrator  as  for  FSK.  Half  the  sum  of  the  inputs 
at  times  t  and  t  -  T  appears  at  time  t  at  one  output  of  the  Mach-Zehnder  filter,  while  half 
the  difference  appeeirs  at  the  other  output.  Thus  depending  on  whether  the  signal  phase  was 
shifted,  the  information  signal  will  appear  only  at  one  of  the  outputs. 

Due  to  the  importance  of  optical  amplifiers  there  have  been  many  recent  of  papers  analyzing 
such  systems,  often  using  approximations,  e.g.  [12],  [20],  [10]  [17],  [8],  [4],  [18],  [6],  as  well  as 
many  other  references  cited  therein. 

If  the  amplifier  spontaneous  noise  dominates  the  receiver  shot  eind  thermal  noises,  the 
performance  of  such  an  optical  system  will  be  identical  to  that  of  a  radio  system  with  square 
law  detection.  This  is  a  well  known  problem  in  communication  theory  [14]  and  the  goal  of  this 
paper  is  to  give  a  self-contained  account  of  the  key  exact  results  in  a  format  appropriate  for 
optical  communication.  We  also  give  some  exact  results  that  can  be  obtained  when  the  shot 
noise  is  considered,  using  the  theory  described  in  [2]. 

2  Analysis 

The  dimensionality  of  the  space  of  finite  energy  signals  with  a  bandwidth  B  and  a  time  spread 
T  is  about  2BT  4-  1  [5],  [19].  For  mathematical  convenience  we  assume  that  2BT  4-  1  is  an 
even  integer,  2M.  For  large  M  one  can  accurately  write  the  filtered  noise  process  at  the 
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photodetector  (in  each  polarization)  over  an  interval  of  time  T  as  ni<f>i{t).  Here  the  <j)is 
are  orthonormal  functions  over  the  pulse  interval.  The  n^’s  are  zero  mean  independent  Gaussian 
random  variables  with  variance  No  1 2.  If  an  information  signal  is  present  at  the  photodetector 
in  the  filter  band  we  can  expand  it  in  the  same  basis,  and  write  it  as  J2i=i  The  energy 

E  of  the  signal  satisfies  E  =  •Sj  >  assuming  that  the  signal  passes  undistorted  through  the 
filter.  This  will  be  the  case  if  M  is  larger  than  1,  as  already  assumed.  The  signal  to  noise 
ratio  E/Nq  can  be  regarded  as  the  number  of  signal  photons  at  the  input  of  the  ideal  high  gain 
optical  amplifier  that  produces  the  noise. 

For  convenience  we  use  units  such  that  the  photodetector  has  unit  gain.  The  integral  of  the 
output  of  the  photodetector  is  then  given  by 


.X  2M  2M 

J  {^{si  +  ni)<i>i{t))^dt  =  J2{ni  +  Si) 

i=l  i=\ 


(1) 


By  completing  the  square  in  the  integral,  it  is  easy  to  see  that  the  Laplace  Transform  of  the 
probability  density  x  is  given  by 


Fe{s)  =  F;[exp(-sa:)] 


1  .  sE  .  V  1 

(1  + 


(2) 


The  subscript  of  the  function  denotes  the  signal  energy.  Fe{-)  is  the  transform  of  a  non-central 
Chi  Square  distribution  with  2M  degrees  of  freedom  [1],  [14].  Note  that  it  does  not  depend 
on  the  individual  s^’s,  but  only  on  the  smn  of  their  squares,  which  is  E.  Tliis  is  due  to  the 
spherical  symmetry  of  the  noise. 

Inverting  Fe{s)  yields  the  probability  density 

1 


-  .  ,  ^  ,  x\  ^  ,  x  +  E,^  ,^\fxE. 

exp(— j^)/M-.(2-j^),  »  >  0 


(3) 


Here  /„  denotes  the  nth  modified  Bessel  function  of  the  first  kind  .  By  differentiating  Fe{s) 
one  finds  that  the  mean  of  the  distribution  is  MNo  -f  E  and  that  its  variance  is  MN^  +  IENq. 
The  two  components  of  the  variance  are  often  denoted  by  “noise/noise  beat”  and  “signal/noise 
beat”  respectively. 
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If  no  signal  is  present  then  E  —  0  and 


{l  +  sNo)^' 

This  is  the  transform  of  the  familiar  central  Chi  Square  density  with  2M  degrees  of  freedom 
(also  called  Gamma  or  Erlang)  [1],  [19],  [14],  with  probability  density 

f  Ij-i  ^  exp(-a:/iVo)  ^ 

These  statistics  will  be  used  in  obtjiining  the  performances  of  dilFerent  modulation  formats. 

When  M  =  1  the  previous  result  are  outside  the  range  of  validity  of  the  model.  However 
they  then  give  the  statistics  of  ideal  direct  detection  systems  with  an  optical  matched  filter, 
and  of  coherent  reception  systems  with  an  IF  matched  filter  and  envelope  detection. 

3  ASK  modulation 

In  the  case  of  ASK  with  equcdly  likely  signals,  the  ON  signal  will  have  energy  2£,  where  £ 
denotes  the  average  signal  energy.  The  receiver  compares  the  output  from  the  integrator  with 
a  threshold  7,  deciding  ON  if  the  output  is  greater  than  7.  The  probability  that  an  ON  is 
decided  when  an  OFF  is  sent  is  given  by 

a  =  /”  M.)d.  =  g  I  (^J  (4) 

Similarly  the  probability  that  an  OFF  is  decided  when  an  ON  is  sent  can  be  expressed  in 
terms  of  the  generalized  Marcum  Q  ftmction  of  order  M, 

Pi  =  f2s{x)dx  =  1  -  Qm{i/4£INo,  y27/iVo) 

[14].  QM{o.,b)  is  defined  as  /a^~^  exp{~{x^  a^)/2)lM-i{o.x)  dx.  Nimierical  methods 

to  compute  it  as  well  as  many  references  to  papers  about  it  appear  in  [16]. 

The  probability  of  error  is  given  by  =  {Pq  +  Pi)l2.  The  optimal  threshold  7  satisfies 
the  equation  /o(7)  =  /2£(7)  or  (2£:7)(^-0/2  =  -  1)!  exp(-2£:/iVo)/M-i(2\/2l^/No). 
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This  can  be  solved  iteratively,  but  it  is  as  easy  to  directly  search  for  the  value  of  7  that 
minimizes  Pg.  Using  the  asymptotic  expansion  of  Im  reveals  that  the  optimal  normalized 
threshold  (7  -  MNo)f2€  approaches  1/4  for  laxge  £ /{NqM^). 

3.1  Gaussian  Approximation 


Because  a:  is  a  sum  of  2M  independent  rmidom  variables  one  is  tempted  to  assume  that  its  prob¬ 
ability  distribution  is  Gaussian.  K  no  signal  is  present  the  mean  of  x  is  MNq  and  the  variance 
MNq.  If  a  signal  is  present  the  mean  is  MNq  +  2£  while  the  variance  is  MN^+ASNo-  Po  and  Pi 
can  be  written  as  Q{{'j  —  MNo)j \JmNq)  and  Q{{MNo  +  2£  —  ‘y)f  +  4t£No)  respectively. 

Here  Q  denotes  the  complementary  distribution  function  of  a  zero  mean,  miit  variance  Gaussian 
random  variable.  The  optimal  threshold  can  be  foimd  nmnerically.  Customarily  [13]  the  thresh¬ 
old  is  set  so  that  Pq  and  Pi  are  equal,  or  7  ==  MNq  +  26 ij MN^  j ( ^ MN^  ij MN^  -j-  4£No). 
We  then  have 


Pe  ^ 


Q{ 


2e/No 


VM+,/M  +  4€/No 


) 


(5) 


3.2  Results 


In  figures  1  and  2  we  show  the  exact  and  approximate  values  of  Pq,  Pi  and  Pg  as  functions  of  the 
normalized  threshold.  We  see  that  they  do  not  match  at  aU,  and  that  the  optimal  thresholds 
are  very  different.  Using  the  threshold  predicted  by  the  Gaussian  approximation  would  result 
in  dismal  actual  performances.  However  the  exact  probability  of  error  at  the  optimal  threshold 
is  imexpectedly  close  to  the  approximate  probability  of  error  at  its  best  threshold. 

The  figiues  also  show  that  the  minimum  of  the  approximate  Pg  occurs  close  to  the  point 
where  Pq  —  Pi,  validating  the  customary  choice  of  threshold  referred  to  above. 

The  exact  Pq  and  Pi  are  also  very  close  to  each  other  when  the  threshold  minimizes  Pg. 
This  is  known  for  M  —  1  [15],  but  many  authors  use  a  Gaussian  approximation  for  v'T  when  a 
signal  is  present  and  incorrectly  claim  that  Pq  dominates  Pi . 

In  figure  3  we  plot  the  exact  and  approximate  values  of  the  probability  of  error  as  functions 
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of  £ INq^  with  M  as  a  parameter.  The  Gaussian  approximation  and  the  exact  method  agree 
within  a  factor  of  5  when  Pe  c:;  10~®,  confirming  the  coincidence  just  noted.  Again  the  normal¬ 
ized  thresholds  differ  widely  as  seen  in  figure  4.  Some  recently  published  approximations  [8] 
underestimate  the  exact  error  probability,  making  the  Gaussian  approximation  appear  worse 
than  it  really  is.  They  also  fail  to  predict  the  optimal  threshold  by  a  wide  margin. 

By  invertmg  the  customary  Gaussian  approximation  (5),  the  signal  to  noise  ratio  necessary 
to  give  a  probability  of  error  is  found  to  be 

+  (6) 

where  q  is  such  that  Pe  =  Qi^)  [4].  This  curve  (for  Pe  =  10~®,  i.e.  g  =  6)  is  displayed  in 
figure  5,  together  with  the  exact  values.  One  sees  that  E/Nq  values  produced  by  the  Gaussian 
approximation  are  too  large  by  about  3  or  4  “photons”,  as  already  noticed  in  [17]  and  [18]. 

The  curve  given  in  [3]  follows  our  exact  answer  for  small  M,  but  it  is  too  low  by  about  5 
“photons”  for  M  =  40.  Curves  such  as  those  in  figure  5  but  expressed  in  dB  with  respect  to 
the  value  for  M  =  1  are  known  as  “noncoherent  combining  loss”  in  the  conmiunication  theory 
literature. 

Using  the  value  of  M  given  by  the  previous  formula  reveals  that  the  noise/noise  beat  MNq 
and  the  signal/noise  beat  term  4£^iVo  become  equal  when  M  =  (2g(l  +  V^))^,  or  M  =  839  for 
g  =  6,  corresponding  to  Pg  =  10“®. 

As  noted  by  [8],  if  the  polarizer  is  omitted  then  twice  as  many  noise  modes  are  admitted  on 
the  photodetector.  Thus  the  probability  of  error  and  the  optimal  threshold  of  such  a  system 
using  a  filter  with  2BT  +  1  =  2M'  can  be  read  on  figures  3,  4  and  5  with  M  —  2M'  as  the 
parameter  value.  This  also  holds  for  the  other  modulation  formats  considered  below. 

4  FSK  and  orthogonal  modulations 

In  the  case  of  orthogonal  modulation  where  each  signal  has  energy  £  one  filter  output  will 
contain  only  noise,  while  the  other  has  both  signal  and  noise.  If  the  integral  of  square  of  the 
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latter  output  is  u,  an  error  occurs  if  the  integral  of  the  square  of  the  output  of  the  first  filter  is 
greater  than  u.  As  in  (4)  this  event  has  probability  P(u)  =  exp(— tt/iVo)  The 

probability  of  error  Fg  is  obtained  by  taking  the  average  of  P(u),  recalling  that  u  has  density 
fs(')  given  in  (3).  One  sees  that  the  results  can  be  expressed  in  terms  of  the  first  M  derivatives 
of  Fs(s)  =  £J[exp(— su)]  given  in  (2),  yielding 

Af-l 


E 

i=0 


(-1)‘ 


d’ 


ds^ 


■Feis/No) 


J  S~l 


One  checks  by  induction  on  i  that 

d^Feis)  _  il{-NoY 


dsi  (i  +  siVo)^+‘  1  +  5iVo  ^  iVo(  1  +  siVo) 
where  Lf(x)  is  the  generalized  Laguerre  polynomial  given  by 

A  / D  +  i](-x)^ 


) 


Tf  (a:)  = 


)  =  E 

k=0  \ 


n  +  kl  ki 


After  replacing  in  (7)  and  interchanging  the  order  of  summation 


1  £ 

^.  =  pr«p(-2]^)E« 


2Nc 


where  c,  =  l/k\  [7],  [14]. 

4.1  Gaussian  Approximation  and  results 


(7) 


(8) 


(9) 


The  variance  of  the  difference  of  the  integrals  is  simply  the  sxim  of  the  individual  variances  (as 
the  noises  are  statistically  independent),  or  2MNq  +  2£No-  The  Gaussian  approximation  thus 
yields  Pg  ^  Q{£/^2MN^  -f  2£No).  Values  of  Pe  given  by  the  exact  and  approximate  formulas 
are  plotted  in  figure  6.  One  sees  that  the  Gaussian  approximation  is  inaccurate  for  low  M,  e.g. 
by  2.6  dB  for  M  =  1  and  Pg  =  10-». 

Using  the  approximation,  the  signal  to  noise  ratio  o/A'o  necessary  to  give  a  probability  of 
error  Pg  =  Q{q)  is 

£/No  +  2M  (10) 
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This  expression  as  well  as  the  exact  values  for  =  10“®  are  displayed  in  figure  5.  One  sees 
that  FSK  is  only  slightly  worse  them  ASK  when  M  =  1  (for  Pg  =  10~®  EjN^  is  40  vs.  38.3  ). 
Comparing  (10)  with  (6)  reveals  that  ASK  has  a  1.5  dB  average  energy  advantage  over  FSK 
when  M  is  very  large.  The  value  of  M  where  the  noise/noise  beat  term  2MNq  is  equal  to  the 
signal/noise  beat  2SNo  is  given  by  M  =  iq^.  When  5  =  6  this  occurs  when  M  is  as  small  as 
144. 

5  DPSK  modulation 

In  the  case  of  DPSK  it  is  necessary  to  expand  the  output  of  the  broadband  filter  over  time 
intervals  (— T,  0)  and  (0,T)  as  both  will  influence  the  output  of  the  Mach-Zehnder  over  time 
(0,T).  We  use  the  same  basis  functions  as  before  and  denote  by  n\  and  s'^  the  coefficients  of 
the  expansion  on  (— T,  0)  and  by  n"  and  s'/  those  on  (0,T). 

Taking  into  account  the  field  scaling  in  the  Mach-Zehnder  filter,  the  integrals  of  the  power 
of  the  two  inputs  wiU  be  respectively  ((n-  +  nf)/2  -f  (s(-  -f  s'/)/2)^  and  “  K')/^  + 

-  ^n/2)^. 

The  correlation  of  and  n'/  is  weak  as  the  correlation  of  the  noise  process  extends  only  over 
time  T/M.  Assuming  it  is  0  ^  all  the  noise  components  {n'^  ±  n(')/2  are  mutually  independent 
and  have  variance  iVo/4,  half  of  what  they  are  in  the  case  of  FSK.  As  in  FSK  one  of  the  signal 
energies  is  0  and  the  other  is  5,  depending  on  the  data  bit.  We  can  thus  conclude  that  DPSK 
has  a  3  dB  advantage  over  FSK  for  all  values  of  M,  ignoring  error  propagations.  This  applies 
to  Gaussian  approximations  as  well. 

The  signal  to  noise  ratio  needed  by  DPSK  to  achieve  Pg  =  10“^  is  displayed  in  figure 
5  as  a  function  of  M.  It  is  not  always  recognized  [11]  that  the  penalty  incurred  by  DPSK 
reception  with  a  broadband  filter  and  an  integrate  and  dump  postdetection  filter  is  related  to 
the  “noncoherent  combining  loss”. 

Our  previous  results  show  that  for  small  M  DPSK  has  about  a  3  dB  advantage  over  ASK 

*I(.  is  0  if  (lie  n'i  and  n’/  are  obtained  l>y  sampling  at.  tiie  Nyi|nist  rale 
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(with  respect  to  average  energy),  but  this  advantage  is  reduced  for  large  M.  Other  comparisons 
with  ASK  and  FSK  appear  in  figure  9  and  are  discussed  in  the  next  section  which  considers 
the  effect  of  the  shot  noise. 

6  Shot  Noise 

In  tliis  section  we  consider  the  fact  that  the  photodiode  is  not  a  perfect  squarer.  We  use 
the  semi-classical  model  where,  conditioned  on  the  envelope  of  the  input  field,  the  photodiode 
generates  electrons  following  a  nonhomogeneous  Poisson  process  with  rate  equal  to  the  square 
of  the  field  envelope.  We  choose  energy  imits  such  that  it  tabes  an  average  of  one  unit  of  energy 
to  generate  a  photoelectron. 

The  total  number  y  of  photoelectrons  generated  over  a  bit  time  by  the  photodiode,  condi¬ 
tioned  on  the  X  defined  in  (1)  has  a  Poisson  distribution  with  mean  x.  Thus  the  conditional 
Z-transform  of  the  distribution  of  y  given  x  is 

E[z^\x]  =  exp((z  —  !)«),  \z\  <  1 


The  expectation  over  x  of  that  last  expression  is  just  the  Laplace  trcinsform  of  the  Clu-Square 
density  of  x  given  in  (2),  evaluated  at  s  =  1  —  z: 


C.(z)  =  =  ^.,(1  -  i.|  < 


Inverting  this  transform  (by  using  (8))  yields  a  Laguerre  distribution  [2] 


9E{k) 


1  d’^Gsiz) 


k]  dz^ 


N, 


E 


E 


z=0 


The  mean  of  y  is  MNq  -f  E,  i.e.  the  mean  of  x,  but  the  variance  MNo{l  -f-  Nq)  -f  E{1  +  2Nq) 
has  increased  by  a  quantity  equal  to  the  mean.  This  increase  corresponds  to  the  shot  noise;  it 
is  relatively  large  when  No  is  small.  If  F  =  0  the  value  of  the  Laguerre  polynomial  in  (7o(-)  is 


distribution  is  negative  binomial. 
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Again  evaluating  the  performances  of  various  modulation  formats  is  straightforward.  For 
ASK  the  receiver  should  compare  y  with  a  threshold  7  to  make  its  decision.  The  optimal  value 
of  the  threshold  is  the  largest  integer  7  that  satisfies  the  inequality  go{‘j)  >  g2s{lf)- 
The  conditional  error  probabilities  are 


and 


Pi  =  J232e{i) 

t=0 


N, 


M-l 


i  +  iVo  ^ 


7  +  i 
.  7  . 


1  +  iVo 


(11) 


The  previous  equality  has  a  simple  justification  using  Bernoulli  trials  with  success  probability 
1/1  +  Nq.  The  first  expression  is  the  probability  that  the  nmnber  of  failures  until  M  successes 
exceeds  7.  The  second  expression  is  the  probability  that  less  than  M  successes  occur  before  the 
7  +  I’st  failure.  The  exact  error  probability  is  Pe  =  {Pq  +  Pi)l2.  The  Gaussian  approximation 
yields  in  this  case 


P  ^n( _ ^ 

V VM(1  +  1/iVo)  +  VM(1  +  1/iVo)  +  4£/Noil  +  l/2iVo) ) 

The  approximate  £/No  necessary  to  achieve  Pe  =  Q{q)  is 

The  same  problem  has  recently  been  considered  in  [6]  leading  to  conclusions  similar  to 
ours  regarding  the  inaccuracy  of  the  Gaussian  approximation  for  M  —  1. 

In  figure  7  we  show  the  necessary  £/No  to  achieve  Pg  —  10~®,  for  various  values  of  No-  To 
interpret  these  curves  recall  that  an  ideal  amplifier  of  gain  G  followed  by  an  ideal  photodiode 
would  result  in  ATq  =  G  —  1  in  our  system  of  units.  Nq  is  reduced  by  losses  between  the  amplifier 
and  the  photodiode  (e.g.  in  the  filters  or  in  the  photodiode).  To  avoid  clutter,  we  have  not 
^Tliat  work  assumes  =  1,  even  though  the  et|ualit.v  only  holds  for  A/  =  1  or  fe  =  0. 
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displayed  the  values  predicted  by  the  Gaussian  approximation.  Their  match  is  slightly  worse 
than  corresponding  values  in  figure  5. 

In  the  case  of  FSK  an  error  certainly  occurs  if  the  photon  count  from  the  “noise  only”  filter 
exceeds  that  from  the  “signal  plus  noise”  filter;  a  tie  occms  when  those  two  counts  are  equal. 
Let  u  express  the  cotmt  from  the  filter  with  the  signal.  Using  (11)  we  express  the  conditional 
probability  of  error  given  u  as 


No  {u  +  l){u  +  2)..{u  +  j)  /  1 

J=0 


1  +  iVo 


j! 


1  +  iVo 


6  denotes  the  delta  function,  its  presence  is  to  take  the  ties  into  account.  The  average  over  u 
of  the  previous  expression  can  be  written  in  terms  of  the  derivatives  of  the  generating  function 
of  u,  Gs(z)  to  yield 


M-l 


J=o 


dPh{x)k{x) 


dx^ 


^  yOeiz)] 

dzJ  I  ^  'h=No/l+No 

^j\  d^~^h{x)  d‘fc(a:) 


t=0 


Of  comse 

and  as  in  (8) 

d^Geiz) 

dzi  ~  {l  +  No{l-z))M+i'-'‘^^l  +  No{l- z)'~^  ^No{l  +  No{l-z)) 
After  replacement  and  interchange  of  summations  one  obtains 


iim 


dx^-'^  dxj 


_exp(^--%-^:,-^?-  YLf  \- 


) 


Pe  =  (fi^)^exp(- 


+  2iVo  ‘ 


l  +  2No 


M-l 
)  cfe 


fe=o 


S{1  +  No) 
No{l  +  2No) 


where 


Cfe 


1 


M-l  / 

S  V(i+iVo)(H-2iVo)^ 


N^ 


'M  +  i-r 

I  i  —  k 


M-i-l 

E 

1=0 


1  +  iVo 


1+1 


i  +  3 


{■^0+2^M-i-l-l) 


For  M  —  1  this  gives  the  known  [2]  result  =  .5exp(— if/(l  +  2No)).  The  previous  formulas 
reduce  to  (9)  when  No  gets  large. 

The  Gaussian  approximation  yields 

e/No  \ 


Pe-Q 


y2A/(l  +  l/No)  +  26 /Noil  +  l/2iVo) 
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The  approximate  CjN^  necessary  to  achieve  Pe  —  Q{q)  is 


eiNo « +  5^)  + 

In  figure  8  we  show  the  necessary  £/Noto  achieve  —  10~®  as  a  function  of  M,  for  various 
values  of  No-  Again  we  have  not  displayed  the  values  predicted  by  the  Gaussian  approximation, 
their  match  is  slightly  worse  than  corresponding  values  in  figure  5. 

The  probability  of  error  for  DPSK  is  the  same  as  for  FSK,  except  that  No  is  reduced  by  a 
factor  of  2.  In  figure  9  we  display  the  ratio  of  the  energy  £  required  by  DPSK  to  that  required 
by  ASK  and  also  to  that  required  by  FSK,  for  Pg  =  10~®.  One  sees  that  DPSK  rapidly  loses 
its  advantage  in  presence  of  shot  noise.  Tliis  is  not  surprising:  in  absence  of  background  noise 
(in  the  “quantum  limited”  regime),  FSK  and  DPSK  have  identical  performances,  but  ASK  is 
twice  as  efficient  in  terms  of  average  energy. 

7  Conclusion 

We  have  given  exact  probability  of  error  expressions  for  optical  signals  in  presence  of  amplifier 
spontaneous  noise  and  photodetector  shot  noise,  and  we  have  compared  these  results  with 
Gaussian  approximations.  The  quality  of  the  approximation  is  questionable  when  it  is  used 
to  evaluate  the  optimal  threshold  for  ASK,  or  the  probability  of  error  for  small  values  of  M. 
Post  detection  receiver  thermal  noise  is  often  also  significant.  Exact  closed  form  results  are  not 
available  for  that  case,  although  receiver  statistics  cein  be  given  in  the  transform  domain.  The 
probability  of  error  can  be  found  by  numerical  inversion,  approximated  by  the  steepest  descent 
method,  as  was  done  recently  in  [9],  or  bounded  by  ChernofF’s  technique,  as  Personick  [12]  did 
two  decades  ago.  One  expects  Gaussian  approximations  to  become  accurate  as  the  thermal 
noise  becomes  more  dominant.  Thermal  noise  will  also  reduce  the  advantage  of  DPSK  over 
ASK  and  FSK. 
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Normalized  Threshold 


Figure  2:  The  exact  and  approximate  probabilities  of  error  Pq,  Pi  and  Pg  as  a  function  of  the 
normalized  threshold  for  M  =  30  and  6 /No  =  65. 
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Figure  3:  The  probability  of  error  for  ASK  as  a  function  of  signal  to  noise  ratio  S/Nq  for 
various  M  values.  For  each  M,  the  exact  result  as  well  as  the  Gaussian  approximation  with 
both  customary  and  optimal  threshold  setting  Eire  shown. 
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Normalized  Threshold 


Figure  4:  The  optimal  normalized  threshold  as  a  function  of  signal  to  noise  ratio  SjNQ  for 
various  M  values.  For  each  M,  the  exact  result  as  well  as  the  Gaussian  approximation  are 
shown,  with  the  approximation  being  generally  too  low. 


Figure  6:  The  probability  of  error  for  FSK  as  a  function  of  signal  to  noise  ratio  f  /iVo  for  various 
values  of  M.  For  each  M,  the  exact  result  as  well  as  the  Gaussian  approximation  are  shown. 
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FigTire  9:  Energy  gain  of  DPSK  vs.  ASK  (solid  curves)  and  of  DPSK  vs.  FSK  (dashed  curves) 
to  insure  Fe  =  10~®  for  various  values  of  as  fimctions  of  M. 
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